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In this review, I discuss briefly how the presence of a cosmological constant in the Universe 
may imply a decoherent evolution of quantum matter in it, and as a consequence a fundamental 
irreversibility of time unrelated in principle to CP properties (Cosmological CPT Violation). In this 
context, I also discuss recently suggested novel possible contributions of massive neutrinos to the 
cosmological constant, which are not due to the standard loop expansion in quantum field theory, but 
rather due to unconventional properties of (some version of) the quantum theory underlying flavour 
mixing. It is also argued that quantum space time foam may be responsible for the neutrino mass 
differences, observed today, and through the above considerations, for the (majority of the) dark 
energy of the Universe in the present era. In the above context, I also present a fit of all the currently 
available neutrino oscillation data, including the LSND "anomalous" experimental results, based on 
such a CPT Violating decoherent neutrino model. The key feature is to use different decoherent 
parameters between neutrinos and antineutrinos, due to the above-mentioned CPT violation. This 
points to the necessity of future experiments, concentrating on the antineutrino sector, in order to 
falsify the model. 



PACS numbers: 



I. INTRODUCTION 



Recent astrophysical observations, using different ex- 
periments and diverse techniques, seem to indicate that 
70% of the Universe energy budget is occupied by "vac- 
uum" energy density of unknown origin, termed Dark 
Energy 0, Q ■ Best fit models give the positive cosmo- 
logical constant Einstein-Friedman Universe as a good 
candidate to explain these observations, although mod- 
els with relaxing to zero vacuum energy (quintessential, 
i.e. involving a scalar field which has not yet reached 
the minimum of its potential) are compatible with the 
current data. 

From a theoretical point of view the two categories of 
Dark Energy models are quite different. If there is a 
relaxing to zero cosmological vacuum energy, depending 
on the details of the relaxation rate, it is possible in gen- 
eral to define asymptotic states and hence a proper Scat- 
tering matrix (S-matrix) for the theory, which can thus 
be quantised canonically. On the other hand, Universes 
with a cosmological constant A > (de Sitter) admit 
no asymptotic states, as a result of the Hubble horizon 
which characterises these models, and hampers the def- 
inition of proper asymptotic state vectors, and hence a 
proper S-matrix. Indeed, de Sitter Universes will expand 
for ever, and eventually their constant vacuum energy 
density component will dominate over matter in such a 
way that the Universe will enter again an exponential (in- 
flationary) phase of (eternal) accelerated expansion, with 
a Hubble horizon of radius Sh oc 1/ VA. It seems that the 
recent astrophysical observations 0, Q seem to indicate 
that the current era of the universe is the beginning of 
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such an accelerated expansion. 

Canonical quantisation of field theories in de Sitter 
space times is still an elusive subject, mostly due to the 
above-mentioned problem of a proper S-matrix defini- 
tion. One suggestion towards the quantisation of such 
systems could be through analogies with open systems 
in quantum mechanics, interacting with an environment. 
The environment in cosmological constant models would 
consist of field modes whose wavelength is shorter than 
the Hubble horizon radius. This splitting was originally 
suggested by Starobinski , in the context of his stochas- 
tic inflationary model, and later on was adopted by sev- 
eral groups yj. Crossing the horizon in either direc- 
tion would constitute interactions with the environment. 
An initially pure quantum state in such Universes/open- 
systems would therefore become eventually mixed, as 
a result of interactions with the environmental modes, 
whose strength will be controlled by the size of the Hub- 
ble horizon, and hence the cosmological constant. In par- 
ticular, for some simple cases, such as conformally cou- 
pled scalar fields |4( in de Sitter spaces it has been shown 
explicitly that the system modes decohere if they have 
wavelengths longer than a critical value, which is of the 
order of the Hubble horizon. Such decoherent evolution 
could explain the classicality of the early (or late, in the 
case of a cosmological constant) Universe phase transi- 
tions p|. The approach is still far from being complete, 
not only due to the technical complications, which force 
the researchers to adopt severe, and often unphysical, 
approximations, but also due to conceptual issues, most 
of which are associated with the back reaction of matter 
onto space time, an issue often ignored in such a context. 
It is my opinion that the latter issue plays an important 
role in the evolution of a quantum Universe, especially 
one with a cosmological constant, and is associated with 
quantum gravity issues. The very origin of the cosmologi- 
cal constant, or in general the dark energy of the vacuum, 
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is certainly a property of quantum gravity. 

Since string theory seems to be the most rigorous and 
most successful approach to quantum gravity, to date, 
encompassing the known quantum field theories of flat 
space times in its low energy limit, I would like to ap- 
proach the problem of the cosmological constant in this 
framework. Critical string theory models, which are 
based on S-matrix theory, at least for their perturba- 
tive formulation, can only accommodate relaxing to zero 
vacuum energy models, which allow for a proper defini- 
tion of asymptotic states, but cannot deal with de Sitter 
Universes [6|. On the other hand, string/membrane the- 
ory models with anti-de-Sitter backgrounds, which admit 
supersymmetry, are consistent. In critical string theory 
an evolution of pure states to mixed do not exist, and 
this is another way of understanding the incompatibility 
of conventional strings with de Sitter Universes. How- 
ever, quantum effects in strings, generated by dilaton tad- 
poles in cr-models formulated in world sheets with higher 
topologies (genera) lead to non-zero contributions to the 
cosmological constant (Fischler-Susskind mechanism) pj. 
Modular divergences in such theories require regularisa- 
tion, which essentially means that small handles in closed 
strings (appropriate for gravity), of size smaller than the 
world-sheet short-distance cutoff, must be integrated out 
in an effective Wilsonian path integral, since such small 
handles could not be distinguished from tree-level world- 
sheet topologies. This corresponds to adding to the tree 
level (T-model action a world-sheet counter-term appro- 
priate for absorbing such modular (small handle) diver- 
gences. The effects of such terms lead to corrections in 
the tree-level graviton world-sheet /3-functions similar to 
those arising in a de Sitter Universe. However, the issue 
of the impossibility of defining a proper S-matrix in de 
Sitter Universe brings up an immediate question on the 
consistency of the approach within string theory. 



II. NON-CRITICAL STRINGS IN DE SITTER 
BACKGROUNDS 

Fortunately there is a way out, which goes beyond the 
above-described Fischler-Susskind mechanism for gen- 
erating cosmological constant contributions. From a 
world-sheet view point, a non-zero cosmological constant 
amounts to contributions to the effective central charge of 
the two-dimensional world-sheet field theory, which thus 
deviates from its (conformal point) critical value. Stringy 
cr-models with a non-zero central charge deficit consti- 
tute the so-called Liouville strings 8] . The path-integral 
world-sheet quantisation of such non-critical strings re- 
quires the introduction of the Liouville mode </>, which is 
an extra world-sheet field whose target-space signature 
(time-like or space-like) depends on the signature of the 
central charge deficit (positive or negative respectively). 
This extra field, plays the role of a new coordinate in tar- 
get space, and its presence is responsible for the restora- 
tion of the conformal invariance of the theory Q . 



An important step towards the physical significance of 
the Liouville mode is the identification Q of its world- 
sheet zero mode <po with the target time, in the super- 
critical (positive central charge deficit) theories. Such an 
identification emerges from dynamics of the target space 
low-energy effective field theory ^(J , in the sense of min- 
imisation of the effective potential. Furthermore it can 
be shown rigorously 9] that under such an identification 
one cannot define a pure-state S-matrix but rather a $ 
-matrix, which is a non-factorisable product of S and S' , 
acting on density matrix mixed states rather than pure 
states. The non factorisability may be attributed to di- 
vergences in the short-distance world-sheet behaviour of 
the cr-model theory. 

With this identification in mind, one may proceed to 
discuss the issue of propagation of quantum matter in a 
de Sitter Background, within such a non-critical string 
framework. By following simple arguments on world- 
sheet renormalisation-group invariance of cr-model quan- 
tities which have target-space physical relevance, it is 
straightforward to arrive at the following master equa- 
tion describing the evolution of string low-energy matter 
in a non-conformal cr-model background: 

p = z[p,ff]+:/3 i; ^,p] : (1) 

where p is the density matrix of string matter exci- 
tations, H is the effective low-energy matter Hamilto- 
nian, g l are background target-space fields, and (3 l are 
their corresponding cr-model renormalisation group (5- 
functions, expressing their scaling under Liouville dress- 
ing :8|. Canonical quantisation for the operators/fields 
g % is possible in Liouville strings 0, as a result of sum- 
ming up higher world-sheet topologies, and thus : • • • : in 
(|T|) denotes appropriate quantum ordering. The quantity 
Qij = 2z 2 z 2 < Vi(z)Vj(0) > is the so-called Zamolod- 
chikov "metric" in the moduli space of the string, a two- 
point correlation function with respect to the vertex op- 
erators Vi corresponding to the deformations of the cr- 
model action from the conformal point S*: 

S a = S*= g* f V % , (2) 

where J s denotes integration over the world-sheet. The 
dot over p in Q denotes differentiation with respect to 
the world-sheet zero mode of the Liouville field, identified 
in this approach with the target time [jj . An important 
note for the compact notation in (J2J is now in order. 
The index i runs over both species of background target 
space fields as well as space-time coordinates. Thus the 
summation over i, j indices in Q corresponds to a sum- 
mation over M, N indices but also a continuous generally 
covariant space-time integration J d d yy/—g, where y de- 
notes a set of d-dimensional space-time coordinates. It 
is important to stress that strings respect general covari- 
ance by construction. For instance, for the case at hand, 
where we are interested in perturbations of the metric 
background gMN, where M,N are target space-time in- 
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dices, one has the correspondence: 
9 l -» 9mn(v) , 

Vi -> V MN (X,y) =d a X M d a X N S^{y-X(a,T)) , 
where a = a, r , 

/ 9% - / d d y^—gg MN {y) f V MN (X,y) , (3) 

For conformal world-sheet backgrounds (5 l = and one 
obtains a normal quantum mechanical equation where 
purity of states is preserved under evolution. When non- 
conformal string backgrounds are present, however, one 
has non-quantum mechanical corrections terms in this 
evolution, which in general may imply decoherence of 
matter, that is evolution of initially pure states to mixed 
ones. In a perturbative derivative expansion (in powers 
of a' , where a' = t 2 s is the Regge slope of the string, 
and £ s the fundamental string length), the lowest order 
graviton (3 function is just the Ricci tensor 

Pmn = a'R M N + ■ ■■ (4) 

where the . . . indicate terms higher order in a', which 
can be ignored in a low-energy (infrared) framework for 
the target space effective field theory, we are interested 
in here. From now on, unless otherwise stated, we shall 
work in units of a' = 1 for brevity. 

Conformal backgrounds in string theory are therefore 
Ricci flat backgrounds in this framework. On the other 
hand, de Sitter backgrounds, for which Rmn oc KguN 7^ 
0, with A > a cosmological constant, obviously vio- 
late this condition, and the excitation of strings in such 
backgrounds can be described, at least perturbatively for 
the (physically relevant) case of small A, so that the de- 
viation from the Ricci flatness is minute, by means of a 
non-critical Liouville dressed a- model. From {Q, J3Jl, the 
evolution equation of low-energy string matter in such a 
background, then, reads in this case 

P = i[p,H} + J 'd d yA: V=99mn [g MN , p] : (5) 

where we took into account the fact that to lowest or- 
der, the Zamolodchikov metric is just the appropriate 
tensorial identity in all sets of indices (M and y). If 
one chooses an antisymmetric ordering prescription, and 
adopts a weak-graviton expansion about flat Minkowski 
space time, gMN — Vmn + h>MN, which seems to be the 
cosmologically relevant case for the present era of the 
Universe, then one arrives at a double commutator struc- 
ture for the decoherence term [Tlj : 

p = i[p,H}+ J d d yk[h MN , [h MN , p}] (6) 

Notice that the decoherence term, which is real, is not 
invariant, due to its structure, under the time reversal 
symmetry t — ► — t (we remind the reader that under 
such a symmetry, the matter Hamiltonian is time rever- 
sal invariant, but the i — » —i. Moreover, since we are 



dealing with small perturbations around flat Minkowski 
space time, the quanta of the gravitational field Hmn can 
be taken to respect the time reversal symmetry). This 
breaking of the time reversal invariance is unrelated in 
principle to properties of matter under the discrete sym- 
metries of Charge (C) and Parity (P), and thus the A- 
induced decoherence term is CPT violating. This is what 
I would call Cosmological CPT Violation , due to the 
global nature of the non-quantum mechanical terms in 
©. Moreover, taking into account that AgMN in de Sit- 
ter spaces may be viewed as a contribution to the stress 
tensor TjjJ^ of the vacuum, one observes that the de- 
coherence term in (jHJl may be considered as a quantum 
version of the (integrated over space time) trace of this 
tensor, thereby being proportional to the global confor- 
mal anomaly of the de Sitter space-time vacuum. 

The above results are in full agreement with the vio- 
lation of CPT in decoherent field theories characterised 
by an evolution of pure to mixed quantum states, which 
we have here due to the presence of the Hubble hori- 
zon Q . Indeed, according to a mathematical theorem by 
R. Wald [12, the CPT operator is not a well-defined 
quantum mechanical operator in field theories where 
there is decoherence, that is evolution of pure to mixed 
states. This leads to a violation of CPT symmetry in 
its strong form, or rather microscopic time irreversibil- 
ity. This may lead to different decoherent parameters 
eventually between particles and antiparticles, reflecting 
the different ways of interaction with the foam between 
the two sectors. We should remark here that, in the case 
of two-state systems, such as two generation neutrino os- 
cillation models, the double commutator terms propor- 
tional to the cosmological constant in may be ex- 
pressed in terms of metric variations (AgMN) 2 between, 
say, neutrino energy eigenstates, expressing back reac- 
tion of neutrino fluctuations onto the space time, as a 
consequence of interaction with the foam |13| . The in- 
duced CPT violation may in general imply, then, metric 
variations of different strength between particle and an- 
tiparticle sectors. 

On the other hand, as stated in it could be pos- 
sible that, despite the strong violation of CPT, a weaker 
form of CPT invariance is maintained phenomenologi- 
cally, in the sense that an observer can always prepare 
pure initial states \<j) >, which could evolve to pure final 
states, >, and for this subset of states the probabilities 
for the transition and its CPT image were equal: 



p(4> -f V) = P(0 -1 V -» H) (7) 



where 9 is the anti-unitary CPT operator acting on pure 
states only. Such an issue can be disentangled experi- 
mentally, and this is the next topic in our discussion. 
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III. CHECKING MICROSCOPIC TIME 
IRREVERSIBILITY IN THE LAB: NEUTRINOS 

The most sensitive, and physically interesting, particle 
probe for quantum-gravity decoherence, to date, appears 
to be the neutrino, for which recently there is mount- 
ing experimental evidence that it carries a non-trivial 
mass. The inequality of neutrino masses among the var- 
ious flavours leads to oscillations, whose properties are 
affected by the above-mentioned decoherent evolution. 

In general, decoherent evolution may be induced by 
other means, such as the presence of ordinary matter, 
which the particle passes through, or the presence of 
quantum space-time foam situation, in which microscopic 
(Planck size) topologically non-trivial metric fluctuations 
may make the ground state of quantum gravity behave as 
a 'medium'. Such contributions are in general indepen- 
dent of the above-described cosmological CPT Violation, 
although, as we shall discuss below, there might be a 
common origin of both quantum space-time foam deco- 
herence and cosmological constant in the following sense: 
according to some speculative scenaria |l4j . a neutrino 
mass difference, and hence flavour mixing, might be the 
result of quantum gravity decoherence, in analogy with 
the celebrated MSW effect [f5j, where contributions to 
the mass difference between neutrino flavours is induced 
as a result of the passage of neutrinos through ordinary 
media. In some approach to the quantisation of flavour 
mixing m field theory [HUll, 

as we shall discuss in the 
next section, one can show that there are non-trivial non- 
perturbative contributions to the vacuum energy of the 
Universe (cosmological constant) from massive neutrinos, 
which are in fact proportional to the (sum) of the mass 
differences, at least in hierarchical neutrino models [l4| . 
According to our discussion above, then, this would im- 
ply cosmological CPT Violation, pointing to the inap- 
plicability of flat-space methods for the quantisation of 
the flavour space, which is thus becoming a curved-space 
time (de Sitter) problem, awaiting solution. 

From a phenomenological viewpoint, one can adopt 
a model-independent approach to arrive at the master 
equation for the time evolution of the neutrino density 
matrix, which could encompass all such foam or cos- 
mological constant CPT Violating effects in a unified 
formalism, without the necessity for a detailed micro- 
scopic knowledge of the underlying physics for the 'envi- 
ronment'. For three generation neutrinos this has been 
done in |13|| . and we next proceed to review briefly the 
situation. 

The mathematical formalism adopted is the so-called 
Lindblad or mathematical semi-groups approach to de- 
coherence [Tsj . which is a very efficient way of study- 
ing open systems in quantum mechanics. The time ir- 
reversibility in the evolution of such semigroups, which 
is linked to decoherence, is inherent in the mathemat- 
ical property of the lack of an inverse in the semi- 
group. This approach has been followed for the study 
of quantum-gravity decoherence in the case of neutral 



kaons in |Igll20j. 

The Lindblad approach to decoherence does not re- 
quire any detailed knowledge of the environment, apart 
from energy conservation, entropy increase and complete 
positivity of the (reduced) density matrix p(t) of the sub- 
system under consideration. The basic evolution equa- 
tion for the (reduced) density matrix of the subsystem in 
the Lindblad approach is linear in p(t) and reads: 

^ = -i[H eS ,p} + ^J2([bnP(t)b]} + {b jP (t),b]}) , (8) 

3 

where H c g is the effective Hamiltonian of the subsystem, 
and the operators bj represent the interaction with the 
environment, and are assumed bounded. Notice that the 
Lindblad part cannot be written as a commutator (of a 
Hamiltonian function) with p. Environmental contribu- 
tions that can be cast in Hamiltonian evolution (commu- 
tator form) are absorbed in H c e. 

It must be noted at this stage that the requirement 
of complete positivity, which essentially pertains to the 
positivity of the map p(t) as the time evolves in the case 
of many particle situations, such as meson factories (two- 
kaon states (^-factory), or two-B-meson states etc.), may 
not be an exact property of quantum gravity, whose in- 
teractions with the environment could be non linear 9]- 
Nevertheless, complete positivity leads to a convenient 
and simple parametrization, and it has been assumed 
so far in many phenomenological analyses of quantum 
gravity decoherence in generic two state systems, such as 
two-flavor neutrino systems pll I22I 1231 ] . 

Formally, the bounded Lindblad operators of an TV- 
level quantum mechanical system can be expanded in a 
basis of matrices satisfying standard commutation rela- 
tions of Lie groups. For a two-level system [ijj, H3 such 
matrices are the SU(2) generators (Pauli matrices) plus 
the 2x2 identity operator, while for a three level sys- 
tem |2 ll. which will be relevant for our purposes in this 
article, the basis comprises of the eight Cell-Mann SU (3) 
matrices Aj , £ = 1, ... 8 plus the 3x3 identity matrix 

3 Let J^, (i = 0, . . . 8(3) be a set of SU(3) (SU(2)) gener- 
ators for a three (two)-level system; then, one may expand 
the various terms in JSJ in terms of to arrive at the 
generic form: 

ij V 

M,f = 0, ...iV 2 -l, i,j = l,...N 2 - 1 (9) 

with N = 3(2) for three(two) level systems, and fijk 
the structure constants of the SU(N) group. The re- 
quirement for entropy increase implies the hermiticity 
of the Lindblad operators bi, as well as the fact that 
the matrix C of the the non-Hamiltonian part of the 
evolution has the properties that Co^, = C^o = 0, 

Aj = jEtAm^m'^ fimkfekj, with the notation bj = 
H,fi bj Jp,- 
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In the two-level case of |l9| the decoherence matrix 
Cfj, v is parametrised by a 4 x 4 matrix, whose non van- 
ishing entries are occupied by the three parameters with 
the dimensions of energy a, f3, 7 with the properties men- 
tioned above. If the requirement of a completely positive 
map p(t) is imposed, then the 4x4 matrix C becomes 
diagonal, with only one non vanishing entry occupied by 
the decoherence parameter 7 > [23 ■ 

In the CPT Violation feature of space-time foam 
decoherence, has been taken into account for the neu- 
trino oscillation case, by assuming different decoherence 
parameters between particle and antiparticle sectors. Be- 
low we shall use the barred notation for the antiparticle 
sector quantities. Notice that this is possible in neu- 
trino oscillations because we are dealing with oscillations 
among flavours separately between particle and antipar- 
ticle sectors, e.g. we shall be interested in probabilities 
P Va ^, Vfj , or P Qoi _,p , where a,f3 are neutrino flavours. In 
contrast, in neutral meson systems one is deal- 

ing with oscillations between particle antiparticle sectors 
(e.g. K° — > K° for Kaons, etc.), and hence the relevant 
decoherent evolution contains only one set of parameters 
(a, /3, 7) for both sectors. 

The extension of the completely positive decoherence 
scenario to the standard three-generation neutrino oscil- 
lations case requires formally the adoption of the three- 
state Lindblad problem. The relativistic neutrino Hamil- 
tonian H c g ~ p 2 + m 2 /2p, with m the neutrino mass, 
has been used as the Hamiltonian of the subsystem 
in the evolution of eq.JSJ). In terms of the generators 
Jfj,, p = 0, ...8 of the SU(3) group, H c g can be ex- 
panded as 0: H oS = ± ^273 (6p 2 + Eti m tj Jo + 
±{Am 2 12 )J 3 + ^ (Am 2 3 + Am 2 23 ) J B , with the obvi- 
ous notation Amfj = mf — m 2 , i,j — 1, 2, 3. 

The analysis of |2J| assumed ad hoc a diagonal form 
for the 9x9 decoherence matrix C in (JSJ: 

[£ M „] = Diag (0, -71, -72, -73, -74, -75, -76, -77, ~7s) 

( 10 ) 

in direct analogy with the two-level case of complete pos- 
itivity [2lll23l |. As we have mentioned already, there is no 
strong physical motivation behind such restricted forms 
of decoherence. This assumption, however, leads to the 
simplest possible decoherence models, and, for our phe- 
nomenological purposes in [l3l | and here, we will assume 
the above form and use it to fit all the available neutrino 
data. It must be clear to the reader though, that such 
a simplification, if proven to be successful (which, as we 
shall argue below, is the case here), just adds more in 
favour of decoherence models, given the restricted num- 
ber of available parameters for the fit in this case. In 
fact, any other non-minimal scenario will have it easier 
to accommodate data because it will have more degrees 
of freedom available for such a purpose. 

In this formalism, the neutrino transition probabilities 



read [Hill: 

PK^)=Tr[p Q (*)/] = 

1 + 1^^*^-1^(0)^ (11) 

i,k,j 

where a, (3 = e, [a,t stand for the three neutrino flavors, 
and Latin indices run over 1, ... 8. The quantities A& are 
the eigenvalues of the matrix M. appearing in the evolu- 
tion @, after taking into account probability conserva- 
tion, which decouples po (t) = -^/2/3, leaving the remain- 
ing equations in the form: dpu/dt = J2j MkjPj- The 
matrices T>ij are the matrices that diagonalise M. [Isj . 
Explicit forms of these matrices, the eigenvalues Afc, and 
consequently the transition probabilities (|11|) . are given 
in H3. 

The important point to stress is that, in generic models 
of oscillation plus decoherence, the eigenvalues Afc depend 
on both the decoherence parameters 7^ and the mass 
differences Amf,-. For instance, Ai = 5 [—(71 +72) — 

\/(72 — 7i) 2 — 4A 2 2 ], with the notation Ay = Am?-/2p, 
i, j = 1,2,3. Note that, to leading order in the (small) 
squared-mass differences, one may replace p by the total 
neutrino energy E, and this will be understood in what 
follows. We note that Afc depend on the quantities $7^ : 



n i2 = 


V (72 


-7i) 2 


-4A? 2 


fiia = 




-74) 2 


-4A 2 3 


^23 = 




-7e) 2 


-4A| 3 



From the above expressions for the eigenvalues Afc, it be- 
comes clear that, when decoherence and oscillations are 
present simultaneously, one should distinguish two cases, 
according to the relative magnitudes of Ay and A7fc; = 
7fc — 7( : v) A^-ij 

I > |A 7 fcf|, and (ii) 2|Ay| < |A-y fc/ |. In 
the former case, the probabilities l|ll|) contain trigono- 
metric (sine and cosine) functions, whilst in the latter 
they exhibit hyperbolic sin and cosine dependence. 

Assuming mixing between the flavors, amounts to ex- 
pressing neutrino flavor eigenstates \v a >, a — e, /i,t in 
terms of mass eigenstates \i>i >, i = 1,2,3 through a 
(unitary) matrix U: \u a >— X)i=i ^ai\ u i This im- 
plies that the density matrix of a flavor state p a can be 
expressed in terms of mass eigenstates as: p a — \v a >< 
v a \ = Ylij UaiUajWi >< Vj\. From this we can deter- 
mine p^ — 2Tr(p Q j7 M ), a quantity needed to calculate the 
transition probabilities (|11|) . 

Due to CPT Violation, as mentioned above, we should 
notice at this stage that, when considering the above 
probabilities in the antineutrino sector, the respective 
decoherence parameters 7^ in general may be different 
from the corresponding ones in the neutrino sector, as 
a result of the strong form of CPT violation. This will 
be crucial for accommodating [u| the LSND result [2f| 
without conflicting with the rest of the available neutrino 
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data. This feature is totally unrelated to mass differences 
between flavors. 

Compatibility of all available neutrino data, including 
CHOOZ and LSND H, can be achieved through a 
set of decoherence parameters in HlOfl such that: all 
the 7^ in the neutrino sector are set to zero, restricting in 
this way all the decoherence effects to the antineutrino 
one where: 

7i = 72 = 74 = 75 
and 

73 = 76 = 77 = 78, (13) 
For the decoherence parameters we have chosen (c.f. 

7l = 2 • 1CT 18 ■ E and 73 = 1 • l(T 24 /£ , (14) 

In the above formulae E is the neutrino energy, and 
barred quantities refer to the antineutrinos. This 
parametrisation guarantees positivity of the relevant 
probabilities. Overall, we have introduced only two new 
parameters, two new degrees of freedom, 77 and 73, 
which, as argued in [13l | was sufficient to account for the 
available experimental data, including the "anomalous" 
LSND results. Furthermore, we have also set the CP vi- 
olating phase of the NMS matrix to zero, so that all the 
mixing matrix elements become real. 

Since the neutrino sector does not suffer from decoher- 
ence, there is no need to include the solar data into the fit. 
We are guaranteed to have an excellent agreement with 
solar data, as long as we keep the relevant mass differ- 
ence and mixing angle within the LMA region, something 
which we shall certainly do. 

As mentioned previously, CPT violation is driven by, 
and restricted to, the decoherence parameters, and hence 
masses and mixing angles are the same in both sectors, 
and selected to be 

Am 2 12 = Am - ; 2 = 7 • 1CT 5 eV 2 , 
Am 2 3 =_Ami 2 = 2.5 • KT 3 eV 2 , 
6>23 = 6*23 = tt/4, 6>12 = 6 12 = .45, 

013 = 013 = -05, 

as indicated by the state of the art analysis. 

At this point it is important to stress that the inclu- 
sion of two new degrees of freedom is not sufficient to 
guarantee that one will indeed be able to account for all 
the experimental observations. We have to keep in mind 
that, in no-decoherence situations, the addition of a ster- 
ile neutrino (which comes along with four new degrees of 
freedom -excluding again the possibility of CP violating 
phases) did not seem to be sufficient for matching all the 
available experimental data, at least in CPT conserving 
situations. 

In order to test our model with these two decoher- 
ence parameters in the antineutrino sector, we have cal- 
culated the zenith angle dependence of the ratio "ob- 
served/ (expected in the no oscillation case)", for muon 
and electron atmospheric neutrinos, for the sub-GeV and 



multi-GeV energy ranges, when mixing is taken into ac- 
count. Since matter effects are important for atmospheric 
neutrinos, we have implemented them through a two- 
shell model, where the density in the mantle (core) is 
taken to be roughly 3.35 (8.44) gr/cm 3 , and the core ra- 
dius is taken to be 2887 km. We should note at this stage 
that a "fake" CPT Violation appears due to matter ef- 
fects, arising from a relative sign difference of the matter 
potential between the respective interactions of neutrinos 
and antineutrinos with ordinary matter. This, however, 
is easily disentangled from our genuine (due to quantum 
gravity) CPT Violation, used here to parametrise our 
model fit to LSND results; indeed, a systematic study of 
such effects [27| has shown that "fake" CPT Violation in- 
creases with the oscillation length, but decreases with the 
neutrino energy, E, vanishing in the limit E — > 00; more- 
over, no independent information regarding such effects 
can be obtained by looking at the antineutrino sector, as 
compared with data from the neutrino sector, due to the 
fact that in the presence of "fake" CPT Violation, but 
in the absence of any genuine CPT breaking, the perti- 
nent CPT probability differences between neutrinos and 
antineutrinos are related, AP9F T = — AP-£T T , where 

ap pa 

APCPT = Paj} _ P __ ; and the Greek ind 

ices denote neu- 
trino flavors. These features are to be contrasted with 
our dominant decoherence effects 7 y 1 (|14|l . proportional 
to the antineutrino energy, E, which are dominant only 
in the antineutrino sector. For the same reason, our ef- 
fects can be disentangled from "fake" decoherence effects 
arising from Gaussian averages of the oscillation proba- 
bility due to, say, uncertainties in the energy of the neu- 
trino beams |28j . which are the same for both neutrinos 
and antineutrinos. We, therefore, claim that the com- 
plex energy dependence in (|14l) , with both L ■ E and L /E 
terms being present in the antineutrino sector, may be a 
characteristic feature of new physics, with the L ■ E terms 
being related to quantum-gravity induced (genuine) CPT 
Violating decoherence. 

The results are shown in Fig. 1 (c), where, for the 
sake of comparison, we have also included the experi- 
mental data. We also present in that figure the pure 
decoherence scenario in the antineutrino sector (a), as 
well as in both sectors (b). For completeness, we also 
present a scenario with neutrino mixing but with deco- 
herence operative in both sectors (d). The conclusion 
is straightforward: pure decoherence is wildly excluded, 
while decoherence plus mixing provides an astonishing 
agreement with experiment. 

As bare eye comparisons can be misleading, we 
have also calculated the x 2 value for each of the 
cases |l]|].From this analysis it becomes clear that the 
mixing plus decoherence scenario in the antineutrino 
sector can easily account for all the available experi- 
mental information, including LSND. It is important to 
stress once more that our sample point was not obtained 
through a scan over all the parameter space, but by an ed- 
ucated guess, and therefore plenty of room is left for im- 
provements. On the other hand, for the mixing-only/no- 
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FIG. 1: Decoherence fits, from top left to bottom right: (a) pure decoherence in antineutrino sector, (b) pure decoherence in 
both sectors, (c) mixing plus decoherence in the antineutrino sector only, (d) mixing plus decoherence in both sectors. The 
dots correspond to SK data. 



decoherence scenario, we have taken the best fit values 
of the state of the art analysis and therefore no signifi- 
cant improvements are expected. At this point a word of 
warning is in order: although superficially it seems that 
scenario (d), decoherence plus mixing in both sectors, 
provides an equally good fit, one should remember that 
including decoherence effects in the neutrino sector can 
have undesirable effects in solar neutrinos, especially due 
to the fact that decoherence effects are weighted by the 
distance travelled by the neutrino, something that may 
lead to seisable (not observed!) effects in the solar case. 

One might wonder then, whether decohering effects, 
which affect the antineutrino sector sufficiently to ac- 
count for the LSND result, have any impact on the solar- 
neutrino related parameters, measured through antineu- 
trinos in the KamLAND experiment |2^. fn order to 
answer this question, it will be sufficient to calculate the 
electron survival probability for KamLAND in our model, 
which turns out to be P Pq ^ P/3 | KamLAND — — -63, in perfect 
agreement with observations. As is well known, Kam- 
LAND is sensitive to a bunch of different reactors with 
distances spanning from 80 to 800 km. However, the 
bulk of the signal comes from just two of those, whose 
distances are 160 and 179 km. These parameters have 
been used to compute the survival probability. It is also 
interesting to notice that in our model, the LSND effect is 
not given by the phase inside the oscillation term ( which 
is proportional to the solar mass difference) but rather by 
the decoherence factor multiplying the oscillation term. 
Therefore the tension between LSND and KARMEN |3fj| 
data is naturally eliminated, because the difference in 
length leads to an exponential suppression. Another po- 
tential source of concern for the present model of decoher- 
ence might be accelerator neutrino experiments, which 
involve high energies and long baselines, and where the 
decoherence L-E scaling can potentially be probed. This, 
however, is not the case. Accelerator experiments typi- 



cally join their neutrino and antineutrino data, with the 
antineutrino statistics being always smaller than the neu- 
trino one. This fact, together with the smaller antineu- 
trino cross section, renders our potential signal consis- 
tent with the background contamination. Even more, 
in order to constrain decoherence effects of the kind we 
are proposing here through accelerator experiments, ex- 
cellent control and knowledge of the beam background 
are mandatory. The new KTeV data |3lJ on kaon decay 
branching ratios, for example, will change the v e back- 
ground enough to make any conclusion on the viability 
of decoherence models useless. After all, the predicted 
signal in our decoherence scenario will be at the level of 
the electron neutrino contamination, and therefore one 
would need to disentangle one from the other. 

Having said that, it is now clear that decoherence mod- 
els (once neutrino mixing is taken into account) are the 
best (and arguably the only) way to explain all the obser- 
vations including the LSND result. This scenario , which 
makes dramatic predictions for the upcoming neutrino 
experiments, expresses a strong observable form of CPT 
violation in the laboratory, and in this sense, our fit gives 
a clear answer to the question asked in the introduction 
as to whether the weak form of CPT invariance is 
violated in Nature. It seems that, in order to account for 
the LSND results, we should invoke such a decoherence- 
induced CPT violation, which however is independent of 
any mass differences between particles and antiparticles. 

This CPT violating pattern, with equal mass spec- 
tra for neutrinos and antineutrinos, will have dramatic 
signatures in future neutrino oscillation experiments. 
The most striking consequence will be seen in Mini- 
BooNE [32, According to our picture, MiniBooNE will 
be able to confirm LSND only when running in the an- 
tineutrino mode and not in the neutrino one, as decoher- 
ence effects live only in the former. Smaller but exper- 
imentally accessible signatures will be seen also in MI- 
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NOS |3^|, by comparing conjugated channels (most no- 
ticeably, the muon survival probability) . We should men- 
tion at this stage that our model is in agreement with the 
strong suppression of dccoherence in the neutrino sector 
expected from astrophysical observations of high energy 
cosmic neutrinos |34| . 

We next remark that fits with dccoherence parameters 
with energy dependences of the form (|14|l imply that the 
exponential factors e Afct in due to decoherence will 
modify the amplitudes of the oscillatory terms due to 
mass differences, and while one term depends on L/E 
the other one is driven by L ■ E, where we have set t = L, 
with L the oscillation length (we are working with natural 
units where c = 1). 

The order of the coefficients of these quantities, 7^ ~ 
10~ 18 ,1CT 24 (GeV) 2 , found in our sample point, im- 
plies that for energies of a few GeV, which are typical 
of the pertinent experiments, such values are not far 
from 7^ ~ Amfj. If our conclusions survive the next 
round of experiments, and therefore if MiniBOONE ex- 
periment |32| confirms previous LSND claims, then this 
may be a significant result. One would be tempted to 
conclude that if the above estimate holds, this would 
probably mean that the neutrino mass differences might 
be due to quantum gravity decoherence. Theoretically 
it is still unknown how the neutrinos acquire a mass, 
or what kind of mass (Majorana or Dirac) they possess. 
There are scenaria in which the mass of neutrino may be 
due to some peculiar backgrounds of string theory for in- 
stance. If the above model turns out to be right we might 
then have, for the first time in low energy physics, an in- 
dication of a direct detection of a quantum gravity effect, 
which disguised itself as an induced decohering neutrino 
mass difference. Notice that in our sample point only an- 
tineutrinos have non-trivial decoherence parameters 77 , 
for i = 1 and 3, while the corresponding quantities in 
the neutrino sector vanish. This implies that there is a 
single cause for mass differences, the decoherence in an- 
tineutrino sector, which is compatible with common mass 
differences in both sectors. If this turns out to be true, it 
could then lead to important conceptual changes in our 
thinking of the problem of particle masses in field theory. 



IV. NEUTRINO MIXING, SPACE-TIME FOAM 
AND COSMOLOGICAL CONSTANT? 

In what follows we will make this assumption, namely 
that decoherence effects, due to interactions with the 
foam, contribute to the Hamiltonian terms in the evo- 
lution of the neutrino density matrix, and result in neu- 
trino mass differences in much the same way as the cel- 
ebrated MSW effect [l5j . responsible for a neutrino mass 
splitting due to interactions with a medium. Indeed, 
when neutrinos travel through matter, the neutral cur- 
rent contribution to this interaction, proportional to - 
Gfit-u/V^, with Gf Fermi's weak interaction constant, 
and n n the neutron density in the medium, is present 



for both v e and (in a two flavour scenario), while the 
charged current contribution, given by \/2GFn e , with 
n e the medium's electronic density, is present only for 
v e . The flavour eigenstates v e ^ can then be expressed 
in terms of fields 2/1,2 with definite masses m^j respec- 
tively with a mixing angle 9, the tilde notation indi- 
cating the effects of matter. The tilded quantities are 
diagonalised with respect to the Hamiltonian of v e ^n hi 
the presence of non-trivial matter media, and one can 
find the following relations between vacuum (untilded) 

and medium parameters'^ sin 2 2(9 ~ sin 2 2fl ( ^ff^i ^j » 

with Am 2 = y/(D - Ato 2 cos26>) 2 + (Am 2 sin26') 2 , D = 
2\f2GFn e k. From this we observe that the medium- 
induced effects in the mass splittings are proportional 
to the electronic density of the medium and in fact, even 
if the neutrinos would have been mass degenerate in vac- 
uum, such a degeneracy would be lifted by a medium. 

To get a qualitative idea of what might happen with 
the foam, one imagines a similar mixing for neutri- 
nos, as a result of their interaction with a quantum- 
gravity decohering foam situation. As a result, there 
are gravitationally-induced effective masses for neutrinos, 
due to flavour dependent interactions of the foam, which 
are in principle allowed in quantum gravity. In anal- 
ogy (but we stress that this is only an analogy) with the 
MSW effect, the gravitationally-induced mass-splitting 
effects are expected now to be proportional to G^nbiik, 
where Gn = 1/Mp is Newton's constant, Mp ~ 10 19 
GeV is the quantum gravity scale, and nbh is a "foam" 
density of appropriate space time defects (such as Planck 
size black holes etc.), whose interaction with the neutri- 
nos discriminates between flavours, in an analogous way 
to the matter effect. Neutrinos, being electrically neutral 
can indeed interact non-trivially with a space time foam, 
and change flavour as a result of such interactions, since 
such processes are allowed by quantum gravity. On the 
other hand, due to electric charge conservation of micro- 
scopic black holes, quarks and charged leptons, cannot 
interact non-trivially with the foam. In this spirit, one 
can imagine a microscopic charged black- hole/anti-black- 
hole pair being created by the foam vacuum. Evapora- 
tion of these black holes (probably at a slower rate than 
their neutral counterparts, due to their near extremal na- 
ture can produce preferentially e + e~ pairs (lighter 
than muons), of which the positrons, say, are absorbed 
into the microscopic event horizons of the evaporating 
charged anti-black hole. This leaves us with a stochasti- 
cally fluctuating (about a mean value) electron (or more 
general charge) density, n^(r), induced by the gravita- 
tional foam, (n£ h (r)) = n j= 0, « h (r)ng h (r')) ^ 0, 
which, in analogy with the electrons of the MSW effect 
in a stochastically fluctuating medium |36j. can interact 
non-trivially only with v e but not with the v^, in contrast 
to neutral black holes which can interact with all types 
of neutrinos (l4j. We assume, of course, that the contri- 
butions to the vacuum energy that may result from such 
emission and absorption processes by the black holes in 
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the foamy vacuum are well within the known limits. For 
instance, one may envisage supersymmetric/superstring 
models of space-time foam, where such contributions may 
be vanishingly small 37] . The mean value (macroscopic) 
part, no, of n^ h (r), assumed time independent, will con- 
tribute to the Hamiltonian part of the evolution of the 
neutrino density matrix, p. In analogy with the (stochas- 
tic) MSW effect |l5ll36|. this part yields space-time foam- 
induced mass-squared splittings for neutrinos: 



(15) 



with non trivial quantum fluctuations (k is the neutrino 
momentum scale). To ensure a constant neutrino mass 
one may consider the case where (n£ h (r)), which ex- 
presses the average number of virtual particles emitted 
from the foam with which the neutrino interacts, is in- 
versely proportional to the (neutrino) momentum. This 
is reasonable, since the faster the neutrino, the less the 
available time to interact with the foam, and hence the 
smaller the number of foam particles it interacts with. 
Such flavour-violating foam effects would also contribute 
to decoherence throu gh t he quantum fluctuations of the 
foam- medium densitv flil l36| , by means of induced non- 
Hamiltonian terms in the density-matrix evolution. Such 
effects assume a double commutator structure [lH l36l l38| 
and are due to both, the fluctuating parts of the foam 
density, as well as the effects of the mixin g ljTo|) on the 
vacuum energy. Indeed, as we discussed in [l4j , and shall 
review briefly below, neutrino flavour mixing may lead to 
a non-trivial contribution to the vacuum energy, in a non- 
perturbative way suggested in Jl71 . Hence, such effects 
are necessarily CPT violating pT|. in the sense of entailing 
an evolution of an initially pure neutrino quantum state 
to a mixed one due to the presence of the Hubble hori- 
zon associated with the non zero cosmological constant, 
which prevents pure asymptotic states from being well 
defined. In that case, CPT is violated in its strong form, 
that is CPT is not a well-defined operator, according to 
the theorem of 01 ■ 

For convenience we shall discuss explicitly the two- 
generation case. The arguments can be extended to 
three generations, at the expense of an increase in 
mathematical complexity, but will not affect qualita- 
tively the conclusions drawn from the two-generation 
case. The arguments are based on the observation 
that in quantum field theory, which by definition re- 
quires an infinite volume limit, in contrast to quan- 
tum mechanical treatment of fixed volume [3^, the neu- 
trino flavour states are orthogonal to the energy eigen- 
states, and moreover they define two inequivalent vacua 
related to each other by a non unitary transforma- 
tion G- x {0,t): |0(% = Gg 1 (t)\0(t)) m: where 6 is 
the mixing angle, t is the time, and the suffix f(m) 
denotes flavour (energy) eigenstates respectively, and 
G^ x {t) 7^ Gg(t) is a non-unitary operator expressed in 
terms of energy-eigenstate neutrino free fields v-i ^ [17| : 

Gg(t) = cxp \ 6 J d 3 x[vl{x)v 2 {x) - v\{x)v\{x)]\ • A rig- 



orous mathematical analysis of this problem has also ap- 
peared in |4jj. As a result of the non unitarity of Gg 1 (t), 
there is a Bogolubov transformation^] connecting the 
creation and annihilation operator coefficients appearing 
in the expansion of the appropriate neutrino fields of the 
energy or flavour eigenstates. Of the two Bogolubov co- 
efficients appearing in the treatment, we shall concen- 
trate on = \V^\e i( - u "' 1+Wk >^ t , with w fcjl = ^k 2 + m 2 , 
the (positive) energy of the neutrino energy eigenstate 
i = 1,2 with mass TOj . This function is related to 
the condensate content of the flavour vacuum, in the 
sense of appearing in the expression of an appropriate 
non-zero number operator of the flavour vacuum 
f (0\a r J.a 
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two-generation scenario 16]. |V^| has the property of 
vanishing for mi — m^, it has a maximum at the mo- 
mentum scale k 2 — mim^, and for k ^> ^m\m,2 it goes 
to zero as: 



\V,\ 



(mi - m 2 ) 5 



k = \k\ > yj m\TTl2 (16) 



The analysis of [l7J argued that the flavour vacuum |0), 
is the correct one to be used in the calculation of the 
average vacuum energy, since otherwise the probability 
is not conserved 41] . The energy-momentum tensor of a 
Dirac fermion field in the Robertson- Walker space-time 
background can be calculated straightforwardly in this 
formalism. The flavour-vacuum average value of its tem- 
poral Too component, which yields the required contribu- 
tion to the vacuum energy due to neutrino mixing, is[l7T|: 



f (0\T 0O \0)f = (p™)%o 

= J2f d3k ^ (/<°l«sHJ°>/ 



8sin 2 / d 3 *(w M + Wfc, 2 )|V^| i 
Jo 



(17) 



where 7700 = 1 in a Robertson- Walker (cosmological) met- 
ric background. The momentum integral in l|17f) is cut-off 
from above at a certain scale, K relevant to the physics of 
neutrino mixing. In conventional approaches, where the 
mass generation of neutrino occurs at the electroweak 
phase transition, this cutoff scale can be put on the elec- 
troweak scale K ~ 100 GeV, but this yields unacceptably 
large contributions to the vacuum energy. An alternative 
scale has been suggested in , namely K ~ v /mim 2 as 
the characteristic scale for the mixing. In this way these 
authors obtained a phenomenologically acceptable value 
for(p™). 

In our case we shall use a different cutoff scale [14j. 
which allows for some analytic estimates of <|17H to be 
derived, as being mathematically consistent with the 
asymptotic form of (|16H , which is valid in a regime of mo- 
menta k 3> ■ s Jm\mi. This cutoff scale is simply given by 
the sum of the two neutrino masses, K = fco = mi + m-z, 
is compatible with our decoherence-induced mass differ- 
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ence scenario, and also allows for a mathematically con- 
sistent analytic estimate of the neutrino-mixing contri- 
bution to the vacuum energy in this framework. For 
hierarchical neutrino models, for which m\ 3> m%, we 
have that fco S> •Jm\mi, and thus, if we assume that the 
modes near the cutoff contribute most to the vacuum en- 
ergy l |l 7(1 . which is clearly supported by the otherwise 
divergent nature of the momentum integration, and take 
into account the asymptotic properties of the function 
Vjj, which are safely valid in this case, we obtain: 

«ac mix > ~ 87rsin 2 0(m x - m 2 ) 2 (mi + ^ 2 ) 2 x 
(V2 + 1 + ) \ oc sin 2 0(Am 2 ) 2 (18) 

in the limit 7712 <C mi. For the (1,2) sector, the corre- 
sponding Am 2 is given by the solar neutrino data and is 
estimated to be Am 2 2 — 10 -5 eV 2 , resulting in a contri- 
bution of the right order. In this way the cosmological 
constant A is elegantly expressed in terms of the small- 
est (infrared, Am 2 ) and the largest (ultraviolet, Mj,) 
Lorentz-invariant mass scales available. It can be ar- 
gued Q that the above choice of the cutoff fco ~ mi+m? 
is consistent with our conjecture on the decoherence ori- 
gin of the neutrino mass difference, due to interaction 
with the foam medium Ijl5|) . Notice that the above way 
of deriving the neutrino-mixing contribution to the dark 
energy is independent of the usual perturbative loop ar- 
guments, and, in this sense, the result (|18f) should be 
considered as exact (non perturbative), if true. 

Some important remarks are now in order. First of 
all, our choice of cutoff scale was such that the result- 
ing contribution to the cosmological constant depends 
on the neutrino mass-squared differences and not on the 
absolute mass, and hence it is independent of any zero- 
point energy, in agreement with energy-driven decoher- 
ence models j23|. For us, it is curved space physics that 
is responsible for lifting the mass degeneracy of neutrino 
mass eigenstates and create the "flavour" problem. This 
is an important point, which may serve as motivation 
(not proof) behind such a cutoff "choice" , which we con- 
jecture is a physical "necessity" . We have argued above 
that such a cutoff "choice" is a natural one from the 
point of view of quantum-gravity decoherence-induced 
mass differences. Detailed models of this fall way be- 
yond the purposes of this brief note. Nevertheless, we 
believe that the above-demonstrated self-consistency of 
this cutoff choice within the remit of our toy model of 
space time foam is intellectually challenging and encour- 
aging for further studies of this important issue. 

The above considerations above were based on the sug- 
gestion of ref. 0| on a Fock-like quantisation of the 
flavour space. There is still controversy in the literature 
regarding the physical meaning of such quantum flavour 
states |4^|, in particular it has been argued that, although 
such states are mathematically elegant and correct con- 
structions, nevertheless they lead to no observable con- 
sequences. However, in view of the results of and of 



the present work, such an argument may not be correct, 
since the mass-squared difference contribution to the cos- 
mological constant is an observable (global) consequence 
of the Fock-like flavour space quantisation. The pres- 
ence of a time independent cosmological constant 118|) 
in the flavour vacuum, which notably is not present if 
one uses instead the mass eigenstate vacuum, implies 
an asymptotic future event horizon for the emerging de- 
Sitter Universe. The flat-space time arguments of 
for the flavour space field theory cannot then be applied, 
at least naively, and the problem of quantisation of the 
Fock-like flavour space is equivalent to the (still elusive) 
quantisation of field theories in (curved) de-Sitter space 
times. In such a case one cannot define properly asymp- 
totic states, and hence a scattering matrix. This will 
lead to decoherence, in the sense of a modified temporal 
evolution for matter states. 

We now remark that in the case of (anti) neutrinos pass- 
ing through stochastic media jsjfl , including space time 
foam [l4|. there are additional contributions to decoher- 
ence, besides the presence of a A-term, which may offer 
a natural explanation of the decoherence parameters of 
|13| . An important source of decoherence in such media 
is due to the uncertainties in the energy E and/or the os- 
cillation length L of the (anti)neutrino beam. In fact, it 
can be shown j2^| that if one averages the standard oscil- 
lation probabilities P Va ^ Vfj over Gaussian distributions 
for E and/or L with a variance c 2 , the result is equiva- 
lent to neutrino decoherence models, in the sense of the 
time dependent profile of the associated probability be- 
ing identical to that of a completely-positive decoherence 
model. One finds for n flavours [271 ]. 

n n 
a=l /3=l,o<6 

(l-cos(2Mm 2 b ) e - 2CT2 ( Am -> 2 ) - 

n n 

2 E E Im ( U LUfiaU ab U* 0b ) X 
a=l b=l,a<b 

sin(2Mm^e- 2ff2 < Am «^ (19) 

where U is the mixing matrix £ = (x), a — 
y/{(x- (x)) 2 = (L/4E)r, and x = L/AE. The result- 
ing form is identical to that of decoherence, as becomes 
evident by noting that the exponential damping factors 
can be written in the form e~ ljL with t = L (c = 1), 
and decoherence parameters of order: 2cr 2 (Am 2 ) 2 = 

jjL, from which jj = Lr^ . There are various 

scenaria that restrict the order of a. In general, the 
acceptable bounds on a may be divided in two ma- 
jor categories, depending on the form of the uncertain- 

ties 0: a, ~ Ax ~ A^ < ^ (ffi + *g*), or 

(Tj < j§j (ijFT? + [fif] 2 ) ■ In three generation 
models the values of the length and energy uncertain- 
ties may vary between flavours, and also between neutri- 
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nos and antineutrinos, as a result of the intrinsic CPT 
violation, hence the subscript j in the above formulae 
(for antiparticle sectors it is understood that j — > j). 
From the above considerations it becomes clear that, for 
L ~ 2E/Am 2 , which is characteristic for oscillations, 
one has decoherence parameters jj ~ (Am 2 j E)r? . It 
is interesting to estimate first the order of decoherence 
induced by conventional physics, for instance decoher- 
ence induced by uncertainties in the measured energy of 
the beam due to experimental limitations. For long base 
line, atmospheric or cosmic neutrino experiments, where 
AL/L is negligible, and AE/E ~ 1 such decoherence pa- 
rameters are found at most of order 7 ~ 10 -24 GeV, for 
the relevant range of energies, and they diminish with 
energy, vanishing formally when E — > 00, which seems 
to be a general feature of conventional matter-induced 
decoherence effects |27j . 

To obtain the decoherence parameters of the best-fit 
model of 01 it suffices to choose for the antineutrino sec- 
tor 7-3 = rg ~ AE/E ~ 1, and r\- = r 2 -~ 10~ 18 -E 2 /Am 2 . 
As seen above, the decoherence parameters exhibiting a 
1/E energy dependence could be attributed to conven- 
tional energy uncertainties occurring in the beam of the 
(anti)ncutrinos. However, the parameters proportional 
to E, if true, may be attributed to exotic physics. The 
fact that rj in general receives contributions from both 
length and energy uncertainties provides a natural expla- 
nation for the different energy dependence of the deco- 
herence parameter of the model of I in the antiparticle 
sector. Indeed, having identified = r-g as decoherence 
induced by 'conventional-looking' energy uncertainties in 
the antineutrino sector, it is natural to assume that the 
71 = 72 oc E decoherence is due to genuine quantum 
gravity effects, increasing with energy, which are associ- 
ated with metric tensor quantum fluctuations. This is 
achieved provided we assume that = r| ~ (AL/L) 2 , 
i.e. these decoherence coefficients are predominantly 
oscillation-length-uncertainty driven, and take into ac- 
count that variations in the invariant length may be 
caused by metric fluctuations, since L 2 = g^L^L", im- 
plying (Ag^) 2 ~ (AL) 2 /L 2 , in order of magnitude. To 
obtain the best fit results of I, then, for L ~ 2E/ Am 2 , 
one needs quantum-gravity induced metric fluctuations 
in the antineutrino sector of order (Ag^) 2 ~ 10~ 18 L • E. 
The increase with energy is not unreasonable, given that 
the higher the energy of the antineutrino the stronger the 
back reaction onto space time, and hence the stronger the 
quantum-gravity induced metric fluctuations. The factor 
10~ 18 may be thought of as being of order E/Mp, with 
Mp ~ 10 19 GeV the Planck mass, although alternative 
interpretations may be valid (see discussion on possible 
cosmological interpretations at the end of the article). 
The increase with L is not uncommon in stochastic mod- 
els of quantum foam, where the decoherence 'medium' 
effects build up with the distance the (anti)particle trav- 
els j^. We also mention at this stage that, apart from 
these effects, in stochastic models of foam there are ad- 
ditional contributions to decoherence, arising from the 



fluctuations of the density of the medium. These too can 
mimic the effects of the best-fit model of in the an- 
tineutrino sector, as discussed in some detail in [l4|. but 
their L-dependence is different from that of the above 
effects. Comparison between short and long baseline ex- 
periments, therefore, may differentiate between the vari- 
ous decoherent contributions. 

At this stage we would like to mention that the above- 
described model of decoherence provides a novel and ex- 
tremely economical mechanism to generate the observed 
baryon asymmetry in the Universe [l4|. through a pro- 
cess of equilibrium electroweak leptogenesis. To this end 
we first recall that, in the analysis of ref. [2(J, deal- 
ing with decoherent evolution in the neutral Kaon case, 
the asymmetries between the semileptonic decays of Kq 
and those of Kq turned out to depend linearly on di- 
mensionless decoherence parameters such as 7 = 7/Ar; 
in the parametrisation of Ellis et al. in [ljj, where 
Ar = Tl — r,5 was a characteristic energy scale asso- 
ciated with energy eigenstates of the kaon system. In 
fact, the dependence was such that the decoherence cor- 
rections to the asymmetry were of order 7 in complete 
positivity scenaria, where only one decoherence parame- 
ter, 7 > was non zero. In similar spirit, in our case of 
lepton-antilepton number asymmetries, one expects the 
corresponding asymmetry to dep end, t o leading order, 
linearly on the quantity 7 = 7/%/ Am 2 , since the quan- 
tity v Am 2 is the characteristic energy scale in the neu- 
trino case, playing a role analogous to Ar in the kaon 
case. The only difference from the kaon case, is that here, 
in contrast to the kaon asymmetry results, there are no 
zcroth order terms, and thus the result of the matter- 
antimatter asymmetry is proportional to the dimension- 
less decoherence parameter 7, which we are going to take 
as the la rger of the two dechorence parameters of our 
model in [13J, discussed in the previous section, namely 
7 — > 71 = 10~ 18 • E/\J Am 2 . In this way, the matter- 
antimatter asymmetry in the Universe is estimated to be 
A = 7^t§ — 7i — 10~ 6 . The numerical coefficient 
10~ 18 on 7 may be thought of as the ratio T/Mp with T 
the temperature, whose value gets frozen at the EW sym- 
metry breaking temperature. Thus, B = n »~ n » ~ 

with n v (n v ) the number density of (anti) neutrinos, rt 7 
the number density of photons and g* the effective num- 
ber of degrees of freedom (at the temperature where the 
asymmetry is developed) which depends on the exact 
matter content of the model but it ranges from 10 2 to 
10 3 in our case. This implies a residual baryon asymme- 
try of order 10 -10 , roughly the desired magnitude. 

Finally, before closing we would like to alert the in- 
terested reader in another possible aspect of the Fock 
space quantisation of the flavour vacuum |17| , advocated 
in 43] . In the case of non-trivial neutrino mixing, the 
flavour Fock states do not satisfy the standard Linear 
Lorentz invariant dispersion relations E 2 = p 2 + m 2 , 
where p is the momentum and m a rest mass, since they 
are superpositions of mass eigenstates satisfying standard 
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dispersion relations but with different masses. The idea 
of [43j is that such states may experience non-linear mod- 
ifications of the Lorentz symme try, of which doubly spe- 
cial relativities is one example |44j . which should guar- 
antee the frame independence of the results. The remark 
we would like to make in connection with this, is that, in 
view of the contributions of the Flavour Fock states to the 
cosmological constant <|18[) . it may be possible to spec- 
ify the non-linear modifications of the Lorentz symmetry 
satisfied by the Fock states, which was not possible for a 
general mixing angle in 43] , by adopting the idea of . 
According to that, the low-energy limit of a quantum 
theory of gravity in a space time with a (positive) cos- 
mological constant A, must be a theory which is invariant 
under a deformed Poincare symmetry, with the pertinent 
(dimensionful) deformation parameter being related to 
the cosmological constant. These arguments are valid as 
long as the theory behaves smoothly in the limit when 
the cosmological constant becomes small as compared to 
the Planck scale (or, in general the scale characteristic 



of the quantum gravity), i.e. in the limit A£p lanck — » 0, 
which is certainly the case of (|18fl . In our case, where 
the neutrino mass differences have been conjectured to 
be the result of a space-time foamy situation, such con- 
siderations become of great importance in determining 
the symmetry structure underlying the non-flat space- 
time quantum field (or string) theory at hand. We hope 
to study such important issues in detail in a future pub- 
lication. 
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